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For lattice QCD with two sea quark flavours we compute 
the static quark antiquark potential V(R) in the regime where 
string breaking is expected. In order to increase statistics, 
we make full use of the lattice information by including all 
lattice vectors R to any given separation R = |R| in the in- 
frared regime. The corresponding paths between the lattice 
points are constructed by means of a generalized Bresenham- 
algorithm as known from computer graphics. As a result, we 
achieve a determination of the unquenched potential in the 
range 0.8 to 1.5 fm with hitherto unknown precision. Fur- 
thermore, we demonstrate some error reducing methods for 
the evaluation of the transition matrix element between two- 
and four-quark states. 



I. INTRODUCTION 

Confinement of quarks is an issue of prime impor- 
tance in the understanding of strong interaction physics. 
While the study of the static quark-antiquark potential 
from simulations of quantum chromodynamics has been 
pushed to rather high accuracy in quenched QCD and 
allows today for a precise determination of the string 
tension, the search for evidence of string breaking from 
Wilson loops in simulations of the full QCD vacuum has 
been futile so far. It seems that the linear rise of the 
static potential continues to prevail even in presence of 
vacuum polarization by quark loops [jl]-^ ■ 

The common explanation for this unexpected finding is 
that present studies cannot really resolve the asymptotic 
time behaviour of the Wilson loops and that a multichan- 
nel analysis including light fermion operators is required 
to achieve sufficient ground state overlap in the available 
time range [Q. In fact, a fully fledged multichannel ap- 
proach has been demonstrated to be a viable technique 
to realize breaking of the string between adjoint sources 
in pure gauge theories or fundamental colour sources in 
Higgs models p-E0[. In the case of full QCD, however, 



it appears overly costly to achieve the required statisti- 
cal accuracy of the generalized Wilson loops which in- 
corporate light quark-antiquark pairs in the initial or 
final states The reason is that one is prevented 

from exploiting the translational symmetry, since this 
would require computation of light propagators P(y, x) 
on any source point location, x (see e.g. Refs. |Ll|]l2l). In 



Refs. Jl^Jl^l stochastic estimator methods with maximal 
variance reduction (so-called all-to-all methods) were ap- 
plied to cope with the fluctuations on the multichannel 
correlator matrix, but failed so far to provide sufficient 
accuracy in the infrared regime. 

On the other hand, if one scrutinizes existing QCD 
potential data from Wilson loops [0-|6) for colour screen- 
ing one will notice that the errors become substantial in 
the region of interest, r ~ 1.2 fm[|]. Thus, there is room 
for suspicion that QCD colour screening has so far es- 
caped detection simply for the lack of adequate precision 
of large Wilson loops. 

The main purpose of the present note is to improve on 
this point by pushing for a high precision 'classical' Wil- 
son loop calculation at large separations in Nf — 2 QCD. 
This is achieved by squeezing maximal information out of 
each vacuum configuration through rotational invariance 
and comprehensive utilization of all possible R-values on 
the lattice. As a result of our "all-i? approach" (ARA) 
we are able to present a long range static potential from 
Wilson loops in unprecedented accuracy. 

In section || we shall describe how we go about in 
the build-up of nonplanar loops to any given R. Section 
III contains our potential analysis from ARA which is 
used on top of existing signal enhancement techniques, 
such as conventional APE smearing |L5| and translational 
averaging. As a first step in the direction of a two channel 
investigation of string breaking we study in section IV 
the noise reduction effect from ARA on the transition 



1 We use capital (lower case) letters for quantities in lattice 
units (physical units). 



1 



correlator between static and static-light quark states, 
QQ and QQqq, respectively. 



II. LOOP CONSTRUCTION 

Our aim is to increase statistics in the regime of colour 
screening, i.e. for large quark antiquark separations, R. 
Obviously, on such a length scale, a given QCD vacuum 
configuration contains plenty of information that can be 
exploited for self averaging and thus for error reduction: 
firstly, one can realize, on a hypercubic lattice, a fairly 
dense set of R values; secondly, for a given large value of 
R = |R|, there are many different three- vector realiza- 
tions R on the lattice. 

We wish to make use of this fact by a systematic con- 
struction of off-axis Wilson loops, W(R, T), in the range 
■Rmin < R < Rmax, with R min = 10 w 1.7 R and 
i? max = 12^3 3.5 Ro, where i?o is the Sommer ra- 
dius (in lattice units) that amounts to tq 0.5 fm [ pH| . 
The construction proceeds by choosing all possible vec- 
tors, R, with integer components C m i n , C m id, and C max 
(in any order of appearance) that obey the inequalities, 



overlap with the physical ground state we would like to 
construct lattice paths that follow as close as possible 
the straight line connecting Rq with Rq. This task can 
be accomplished by a procedure which is known as the 
Bresenham algorithm pl[ | in computer graphics. There 
one wishes to map a straight continuous line between two 
points, say and C = (C max ,C m i n ), onto discrete pix- 
els on a 2-d screen. Then one has to find the explicit 
sequence of pixel hoppings in max- and min-directions 
such that the resulting pixel set mimicks best the contin- 
uum geodesic between points and C. The Bresenham 
prescription is simply to move in max-direction unless a 
step in min-direction brings you closer to the geodesic 
from to C, where we assume | C m ax | 2^ | Cmin | ■ R can 
easily be embodied into a fast algorithm based on local 
decision making only (see Fig. pi), by means of a char- 
acteristic lattice function \ that incorporates the aspect 

ratio Cmax/Cmin- 



-Rmin < -R — C m in + C 



2 

mid 



max — max 



(1) 



where |C min | < |C mid | < |C max |. 

Subsequently, the set of solutions to the constraint 
eq. (Q) is sorted according to the correponding values 
of R. In Table | we display the large number of possi- 
ble R values and vectors, R, obtained in this way, with 
the additional restriction, |C max | < 12. So far, only 
R vectors with (|C max |, |C mid |, |C min |) being multiples of 
(1, 0, 0), (1, 1, 0), (1, 1, 1), (2, 1, 0), (2, 1, 1) or (2, 2, 1) have 
been considered p~7| p0|| . Within the investigated regime, 
10 < R < 12V3, we achieve a gain factor of more than 
eight in terms of the spatial resolution in R (from 21 to 
175 different values). Moreover, the number of different 
R vectors that yield the same distance, R, is increased 
by an average gain factor of more than four! Accord- 
ingly, we find the ARA to reduce the statistical errors on 
potential values by factors around two. 

We shall briefly discuss the construction of the gauge 
transporters connecting the quark and antiquark loca- 
tions, Rq and Rq = Rq + R, that appear within the 
ARA nonplanar Wilson loops. In order to achieve a large 

TABLE I. The number of different solutions to eq. (|l| ob- 
tained by ARA for the interval 10 < R < 12^, \Ri\ < 12. 





#i? values 


#R vectors 


#R vectors/_R-entry 


standard 


21 


302 


14.4 


ARA 


175 


11486 


65.6 



FIG. 1. Illustration of a path construction by the Bresen- 
ham algorithm, for C = (5,3). 

The algorithm in two dimensions looks like this: 



cmax2 := 2*cmax 
cmin2 :— 2*cmin 
chi :— cmin2 -cmax 
FOR i := 1 TO cmax DO 
step in max-direction 
IF chi > THEN 

chi :— chi - cmax2 
step in min-direction 
ENDIF 

chi :— chi + cmin2 
ENDDO 

The generalisation to three dimensions is achieved by 
combining two of these 2d-algorithms with different x's 
for max-mid and max-min in just one loop over the max- 
direction. 

In order to convey an idea about the statistics gain 
inherent in such a systematic approach we have listed 
the number of R vectors constructed in this manner in 
Table |. Note that for plane or space diagonal separa- 
tions, we average over the two or six equivalent paths, 
respectively. 
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III. THE STATIC POTENTIAL AT LARGE R 

We base our analysis on 184 vacuum configurations, 
separated by one autocorrelation length, on L\L T = 
24 3 x 40 lattices of N f = 2 QCD at /3 = 5.6 and 
n sea = .1575 (corresponding to m w /m p = .704(5)), pro- 
duced by the T^L-collaboration on an APE 100 tower at 
INFN. These parameter values correspond to the biggest 
lattice volume at our disposal, L a a « 2 fm. In order to 
minimize finite-size effects and possible violations of rota- 
tional symmetry on the L a — 24 torus, C max has been re- 
stricted to |C max | < 12. The lattice constant a was deter- 
mined from the Sommer radius r = R a a w 0.5 fm ]16| ], as 
obtained in our previous investigation |2(i|l , i?o = 5.89(3). 

Before entering the Wilson loop analysis the configu- 
rations are smoothened by spatial APE- or link smear- 
ing @, 



link — > a x link + staples 



(2) 



followed up by a projection back into the gauge 
group [ pT[ , with 26 such iterative replacements and the 
parameter value, a = 2.3. 

The potential values have been obtained by means of 
single and double exponential fits to smeared ARA Wil- 
son loop data W(R,T) within the range, T min < T < 8. 
We shall quote statistical errors that are obtained by 
jackknifing. At the large R values that are of interest 
in view of screening effects, the quality of the statistical 
signal did not allow to include T values larger than 8. 




FIG. 2. The estimate on the static potential as obtained 
from Wilson loops at k = 0.1575 and (3 = 5.6 on 24 3 x 40 lat- 
tices, extracted from a single exponential fit to Wilson loops 
for 4 < T < 8. 

In Fig. H we display estimates on the potential in the 
range 1.7 < r/ro < 3.5 obtained from a single exponen- 
tial fit with T m j n = 4 (they agree with results from a 
double exponential fit with T min = 1). Note that for 
r > 2.4 tq the data are to be interpreted as strict upper 
limits on the potential. The slope is in agreement with 



the string tension, K = aa 2 = 0.0372(8) = 1.139(4)i?o 2 , 
as quoted in Rcf. from a fit to data obtained at 
smaller r < 2.04 rn. Around the separation r c « 2.3 rn 
the potential energy crosses the expected threshold for 
string breaking, 2mps a = 1.256(13), which is indicated 
by a horizontal error band. The errors quoted are statis- 
tical and are well below 1.5 % for r/ro < 3. We find the 
data to follow perfectly a straight line: a flattening of the 
Wilson loop potential is not visible within the accessible 
T range and present statistical errors. 

To complement this result, we computed the overlaps 
with the ground and first excited states by means of two- 
exponential fits, as displayed in Fig. ||. Again our two 
data sets exhibit linear dependencies on r, with nearly 
opposite slopes such that their sum turns out to be close 
to one. The remainder does only slightly depend on r 
and is of order 10 % . 

We conclude that Wilson loop operators are definitely 
not very well suited to uncover string breaking. To 
achieve the required overlap with four-quark states one 
has to introduce them explicitely into the calculation, in 
form of a coupled channel analysis. As a first step in such 
a program, we shall explore in the remaining part of this 
letter a number of signal enhancement techniques based 
on ARA, in application to the transition matrix element 
between two- and four-quark states. 
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FIG. 3. The ground and first excited state overlaps 
(top and lower data sets, respectively) plotted versus the 
quark-antiquark separation, as obtained from Wilson loops. 
The data result from two-exponential fits in the region 
1 < T < 8. 



IV. NOISE ON THE TRANSITION OPERATOR 

In a two-channel approach one extends the Wilson loop 
vacuum expectation value, Cn, to a 2 x 2 correlation 
matrix C as pictogrammed in Fig. ||. The task is then to 
solve the generalized eigenvalue problem 



C{t + T)VL % = Ai(T)C(*)Ui 



(3) 
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where the eigenvalues connect to the two energy levels 
Ei at large enough t: 



Ai(r) = exp(-£Vr) 



(4) 



Unless C happens to be diagonal, the physical states Uo 
and Ui are mixtures of two- and four-quark states. It is 
obvious that the energy levels are insensitive to changes 
in the normalization of the wavefunctions, U;. Thus, in 
addition to spatial APE link smearing, source smearing 
techniques as known from spectrum calculations are ap- 
plicable (see Fig. ||) and will be tested. 



C = 



( 



\ 
) 



FIG. 4. The two channel correlation matrix C. Wiggly 
lines correspond to light quark propagators, solid lines to 
products of gauge links. 

While the noise level on C\\ could be greatly reduced 
by applying ARA on top of standard volume self averag- 
ing (VSA), the quark propagators P(y,x) (wiggly lines) 
entering the remaining components of the correlation ma- 
trix C prevent the direct use of VSA that requires ma- 
trix inversions on all possible source points, x. Pennanen 
and Michael have tested noisy estimator techniques on 
P for curing this problem but did not succeed to reach 
the accuracy necessary to solve eq. (||) Q. This moti- 
vates us to explore noise controlling strategies based on 
ARA rather than VSA on the transition matrix element 
C2i(R,T). One should keep in mind that ARA can be 
put to work at little extra cost since one inversion on a 
single source x renders P(y,x) for all sink positions y. 



FIG. 5. The source smeared transition matrix element C21. 

In addition to APE link smearing, we have explored 
source smearing as illustrated by the diagram, Fig. |^. 
This is applied to the light propagator source, Xxi an d 
consists of iterative replacement (as in previous heavy- 
light spectrum analyses [2^j24| ), 



Xx -> Xx + a 



±3 

V Y Par '- U 
i=±l 



(5) 



which is repeated 50 times, with weight factor, a = 4.0. 
The source is put to zero subsequently outside the volume 



|r s — x| < R q a, with R q — 5. Note that quark propa- 
gators are computed without link smearing throughout 
this work. 

Fig. |^ shows the i?-dependence of the relative errors 
on C2i(R,T = 1) throughout the stringbreaking region 
and illustrates the performance of various signal enhance- 
ment tools that are added on top of ARA: (a) the circles 
correspond to no source or link smearing; (b) crosses re- 
fer to local sources and smeared links; (c) stars refer to 
smeared quark sources and smeared gauge sinks. 
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FIG. 6. Relative errors on C 2 i(R,T = 1) vs. R with ARA. 
String breaking is expected around R — 13. Open circles: no 
link or source smearing; upright crosses: link smearing only; 
stars: smeared source and smeared gauge links. 

Obviously link smearing helps a lot but leaves us still 
with errors of order 50 % in the region R < 18. It turns 
out that smearing is capable of cutting down noise am- 
plitudes further. We find an additional reduction of more 
than a factor two, to a 20 % level at T = 1. 

Unfortunately, such accuracy does not suffice, as it 
cannot be sustained at larger T values, where one wishes 
to analyse C at the end of the day. Given the dense set 
of R- values available from ARA and in view of the fact 
that B{R) = ln[C2i (R,T)] is a smooth function of R, 
there is opportunity to further improve by filtering the 
sequence, {B(Ri)} — * {Bf(Ri)}. We have chosen a fil- 
ter that weighs each individual jackknife sample with the 
fluctuation calculated on the entire data set, 

B f (R i ) = Nr 1 °7 2B ( R j)> ( 6 ) 

|flj-fl«|<i?/ 



with normalization 



\Rj-Ri\<R f 



We found best results with filter radius Rf = 0.5. Let 
us now consider the transition matrix element at T = 5. 
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The effect of filtering is illustrated in Fig. fj] where we con- 
front the signals from ARA plus link and source smear- 
ings, {B(Ri)}, (upper sequence with large error bars) 
with the filtered one, {Bf(Ri)}, (lower data sequence). 
In order to display the systematics of windowing, we have 
refrained from thinning the latter sequence as one should 
in actual applications. We find a striking noise reduc- 
tion through filtering within the large distance regime, 
10 < R < 18: at T = 3 and T = 5 we encounter errors of 
less than 10 % and 20 %, respectively. 
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FIG. 7. The signal for \nC 2 i at T = 5. The filtered data 
have been vertically shifted by —8 to avoid cluttering. String 
breaking is expected around R — 13. 



V. SUMMARY AND CONCLUSIONS 

Exploiting the dense set of R values available at large 
quark-antiquark separations on the lattice we succeeded 
in improving the precision of the QQ potential from Wil- 
son loops in the string breaking regime by a factor of 
two with respect to standard methods. This enables 
us to analyse Wilson loop data well beyond the point 
where string breaking is expected and corroborates pre- 
vious conjectures that Wilson loops do not bear enough 
overlap with QQqq states for uncovering string breaking 
within the T-range available at present. 

The success of our all-i? approach to the Wilson loops 
in the string breaking regime has encouraged us to carry 
out a feasibility study on error control of the transition 
correlator, C21. By additional use of source smearing 
and filtering techniques, we arrive at reasonable signal- 
to-noise ratios for lnC2i. 

For the final chord in a full two-channel analysis of QQ 
and QQqq states one would have to consider the corre- 
lator C22 which contains both, a connected and a dis- 
connected contribution: the former can be tackled with 
our present techniques and we expect sufficient accuracy 
with a factor of T max more effort than for C21; the latter 
relates to the situation of BB pairs at large separations. 



From a previous BB study |2j| one would anticipate that 
the {BB) correlator is dominated in our i?-range and for 
our purposes by (B) (B) . This is currently being investi- 
gated. 
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